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1 INTRODUCTION 



M83 (NGC 5236) is a nearby (4.5 Mpc; Ijhim et al.ll2003l) 
grand design spiral galaxy SAB(s)c (|de V aucoule urs et al.l 
Il99ll ) in a group containing several bright galaxies including 
Centaurus A (NGC 5128) and NGC 4945. Within the cen- 
tral 20" lies the photometric peak (hereafter referred to as 
the nucleus) offset from the centre of symmetr y of the outer 
bulge isophotes and the global gas kinematics (IWolstencroft 



bulge isophotes and the global gas kinematics (I Wolstencroti 
19881: ICallais et al.lll99ll : lTh"atte et al.ll2000l : ISakamoto eTal 



20041 ') and surrounded by a semicircular starburst region, 
some 3" - 8" from the nucleus. M83 has a close dynamical 
companion, NGC 5253 ('Rogsta d et al.|[l974h which also con- 
tains recent star formation, albeit much more compact than 
that of M 83. The closest passa ge of NGC 5253 occurred 1-2 
Gyrs ago (|Rogstad et al.lll974l ) and so the activity in these 
galaxies is likely driven by internal inflow of gas via the bar 
rather than from direct interaction as the clusters are mostly 
1-10 Myrs in age (jHarris et al.ll200ll . hereafter HOI). 



1.1 Central star formation 



ABSTRACT 

We combine VLT/ISAAC NIR spectroscopy with archival HST/WFPC2 and 
HST/NICMOS imaging to study the central 20"x20" of M83. Our NIR indices for 
clusters in the circumnuclear star-burst region are inconsistent with simple instanta- 
neous burst models. However, models of a single burst dispersed over a duration of 
6 Myrs fit the data well and provide the clearest evidence yet of an age gradient along 
the star forming arc, with the youngest clusters nearest the north-east dust lane. The 
long slit kinematics show no evidence to support previous claims of a second hid- 
den mass concentration, although we do observe changes in molecular gas velocity 
consistent with the presence of a shock at the edge of the dust lane. 
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anomalous young cluster at this old end). They were not 
alone in this conclusion: iHeap et al.1 (Il993l ). using WFPC 
photometry and lUE UV spectroscopy also suggested that 
the age gradient went from south-east to north-west, the 
youngest cl usters being in the south-east end of the arc. 
iHeap et al.l also identified t hat the old s t ellar nucleus was 
not a site of Ha emission. iPuxlev et al.l (|l997l ) performed 
low resolution NIR spectroscopy along a single slit bisecting 
the starburst arc and found evidence for a radial age gradi- 
ent. lElmegreen et al' 1 l|l998l ) found a double ring structure 
in their (J-K) colour images of the centre of M83 and asso- 
ciated each ring with an inner Lindblad resonances (ILR). 



iGallais et al.l (|l99ll ) first presented NIR images of sub- 
arcsecond resolution for the central region of M83 which 
lead to an improved understanding of the morphology: the 
region of young stars is composed of a bright circum-nuclear 
arc o f star clusters , located ~7" from the old nucleus. How- 
ever, [Gallair^er^l] also attempted to date the young clusters 
and (contrary to more recent work, discussed later) sug- 
gested that the clusters were younger on the south-east end 
of the arc and older on the north-west end (but with one 
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However, when HOI made a comprehensive study of the 
star clusters using HST/WFPC2 narrow- and broadband 
data coupl ed with the latest Sta rburst99 population synthe- 
sis models (|Leitherer et al ] | 19991 . hereafter SB99), they iden- 
tified that the tangential age gradient along the starburst arc 
was in fact from north-west to south-east with the youngest 
clusters in the north-west (in direct con fiict with earlier wor k 
but resolving the anomalous cluster in lGallais et al.l|l99ll ). 
HOI also found the majority of the clusters in the arc to have 
ages in the range 5-7 Myrs and found evidence of age gradi- 
ents perpendicular to the starburst arc with younger clusters 
along the periphery; they concluded such gradients were in- 
dicative of an inside-out propagation of star formation. Sub- 
sequently, iBresolin fc Kennicutti (,2002 ) found Wolf-Rayet 
features in UV spectra of clusters near the north-west end 
of the arc and dated this star formation at 5 Myrs (and 
the other regions at 7 Myrs), in rough agreement with the 
estimates in HOI for the similar regions. 
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1.2 Nuclear Rings 

Considering the evidence thus far presented, one would most 
likely conclude that the circumnuclear arc seen at the centre 
of M83 is a nuclear ring. Nuclear rin gs are oft en found at the 
centres of barred spiral galaxies (iButalll986t) as are nu clear 
spirals l|Martini fc Poggeill999l : IPogge fc Martinill2002l) and 
their f ormation is thought to be very similar. lAthanassoulal 
j 19921 ) first showed how the orbits in barred galactic poten- 
tials lead to gas shocking on the leading edges to form dust 
lanes with the shocked gas then flowing along the potential 
to the centre. These shocks have been observed across dust 
lanes, both directly as velocity jumps in gas kinematics and 
indirectl y from increased rad io emission (both discussed in 
detail bv I Athanassoulall 19921 ) . 

However, lAthanassoulal 's early simulations didn't have 
the resolution to properly probe the central regions to de- 
termine exac tly what hap pened to the inflowing gas; this 
w as done by Finer et al.l l|l995i) using the similar models 
to lAthanassoulal but with cyli ndrical coord inates for better 
spatial resolution at the centre. [Piner et ahl found that if the 
barred potential permitted two ILRs, the gas flow builds up 
between them to form adense nu clear ring. This was also 
confirmed by ICombed l|l996l) and iButa fc CombesI (Il996l ) 
who went further to say that even in a potential that ad- 
mits one ILR, the inflowing gas will build up there to form 
a ring. Given the rings and ILRs found by Elmcgrc en et al.l . 
it appeared likely that the arc of young stars in M83 was 
formed from gas building up in this way and leading to star 
formation. 

This ILR picture of nuclear ring formation has per- 
sisted despite re c ent w ork suggesting subtle differences. 
iRegan fc TeubenI ()2003l ) pe rformed an exhaustive number 
of simulations like those of IPiner et al.l for various barred 
potentials and they conclude that it is not the presence of 
ILRs that dictate the formation of a nuclear ring, but the 
existence and properties of an orbit family, namely the X2 
family: if this orbit family exists, so can a nuclear ring and 
the authors find an excellent correlation between the X2 or- 
bit of largest extent along the bar major axis and the ra- 
dius of the nuclear ring. If the X2 family does not exist, 
the authors never find a nuclear ring. There is some dis- 
cussion as to definitions: an ILR is normally defined as the 
locus at which the pattern frequency of the bar equals the 
circular orbital frequency less half the epicyclic orbital fre- 
quency, e.g. 0.^ = 0, — k/2, and in weak bars this same 
definition tr aces the largest ex tent of the X2 orbit family. 
However, as Regan fc TeubenI are careful to point out, in 
strong bars with a large quadrupole moment and low cen- 
tral mass concentrations, this is no longer the case and there 
can be regions whe re Q.i, < Q, — k/2 but no X2 orbits exist. 
iMacieiewskil (|2004l ) has also performed high resolution hy- 
drodynamical simulations of the central regions of barred 
galaxies and shows how nuclear rings and nuclear spirals 
are formed b y the same hydrodynamical processes. Although 
IMacieiewskil continues the old trend in matching ILRs with 
the nuclear spiral and ring phenomenon, he is careful to de- 
fine the ILR as the largest extent of the X2 orbit family, 
implying that the X2 family exi sts when referring to an ILR. 
We will retu r n to t he results in lRegan fc TeubenI (|2003l ) and 
IMacieiewskil l|2004l ). particularly to their predictions of the 
phase space distribution of the gas (velocity and density) in 



the ring and to the trajectories of stars on X2 orbits, formed 
in the ring. 



1.3 Nuclear Disks 

There is clearly overwhelming evidence that gas can ac- 
cumulate in the central regions of barred potentials, al- 
though the present high resolution hydrodynamical simu- 
lations do not predict if and how the gas may give rise 
to star formation or how the outflows from the young 
stars and supernovae would aff'ect the gas flow on the ring. 
Lower resolution N-body simulations including stars, gas 
and star formation and then combined hydrodynamical and 
N-body simula t ions w ith st ar formation were performed by 
IWoziiiak et al. I (|2003| ) and lWozniak fc Champaver^ (|2006lf l 
after Emsellem et al. (|200ll ) observed drops in the velocity 
dispersion at the centr es of double-barred spirals (also seen 
in barred spirals by Ch ung fc Bureaij||2004 ) and couldn't re- 
produce them in simple dynamical models. The simulations 
with gas flows, stellar motions and star formation showed 
that when gas shocks in a barred potential and flows to the 
centre, it can produce a central disk. This disk of cold gas 
then produces young stars which dominate the luminosity 
of the central region and give rise to a very low velocity dis- 
persion because of the low random motion of the stars in 
the disk. An important difference between these simulations 
and the pure hydrodynamical ones discussed earlier is that 
the gas builds up in a disk and not a ring or spiral; subse- 
quent star formation is similarly distributed in a disk and 
not a ring. Whether this difference is a resolution problem 
or something more subtle remains to be seen, although if the 
surface density of the gas was higher in the central region of 
the disk, one would expect star formation to initiate there 
and propagate o utwards, potentially l eaving an outer ring of 
young stars as in lKennev et al.l (Il993l ) and a radial age gra- 
dient in the star formation (co nsistent with the observations 
of M83 made bv lPuxlev et al.l . although HOI report younger 
star clusters around the periphery of the arc, both on the 
inner and outer edges). 



1.4 Polar rings and disks 

However, there are certain complications with the circum- 
nuclear ring interpretation for the centre of M83. Most pre- 
vious studies dis cuss an arc of young stars, not a complete 
ring, apart from ICallais et al.l ( 199ll) who claim to observe 
a complete ring or ellipse of star clusters from thei r sub- 
arcsecond H-band image, and lElmegreen et al.l (fl998) who 
found circum-nuclear rings in their (J-K) NIR images. The 
fact that complete ring structures are seen in the NIR and 
not in the optical is highly suggestive of an obscuration of 
the complete ring by dust; an idea supported by extinction 
estimates we present later in JSl And there may be a sim- 
ple explanation for the obscuration o f one h alf of a ring and 
not the other: ISofue fc Wakamatsul l| 19941 ) argue that the 
north-east dust lane is observed to warp out of the plane 
of the main disk as it approaches the nucleus and they go 
so far as to call it polar; a similar warp of the south-west 
dust lane on the other side of the galaxy might explain why 
it is also not visibl e. The hydro d ynam i cal simulations of 
lAthanassoulal (|l992l ): iPiner et all |l995l ): iRegan fc TeubenI 
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(|2003l l; iMacieiewskil l|2004 ) are all in 2D and are unable 
to permit such warping out of the disk plane. We there- 
fore have no theoretical predictions of gas inflow produc- 
ing warped disks, spirals and rings, although we can see 
similar b ehaviour in other galaxi es such as NGC 5383 (as 
noted in ISandaee fc Bedkg|l994h . If similar hydrodynam- 
ical processes occur in 3D as they do in 2D but the gas 
inflow can build up in a polar or inclined ring out of the 
plain of the galactic disk, we would certainly expect the star 
formation to do likewise and given a highly obscuring and 
dusty galactic disk, for the inclination of M83 to our line- 
of-sight we would only see the near side of the inner warped 
disk and an incomple te ring of young stars, as explain ed by 
ISofue fc Wakamatsul . However, ISofue fc WakamatsrJ never 
comment on the plausibility of the galaxy's inclination giv- 
ing rise to the obscuration of the dark south-west dust lane 
by the luminous bulge. 



1.5 Interloping hidden mass concentrations 

After studying the c entre of M83 with n ear infra-red (NIR) 
long-slit kinematics, iThatte et al.l l|2000l . hereafter TOO) ini- 
tially suggested the possibility of a second nucleus or mass 
concentration after the stellar kinematics across the princi- 
pal nucleus (Slit E in Fig. [31 c utting the arc of yo ung stars 
and the 'bridge' identified in ICallais et al.l Il99lh revealed 
two peaks in the velocity dispersion — one associated with 
the nucleus and another which was 2'.' 7 SW of the nucleus. 
Assuming that the stars are dynamically relaxed in the grav- 
itational potential, the secondary peak in the dispersion is 
best explained by invoking the presence of a second dynam- 
ically h ot but obscured m ass concentration (TOO). Subse- 
quently, iMast et al.l ()2006l ) performed optical integral field 
spectroscopy (IFS) on the central region and linked a strong 
velocity gradient in the Ifo; velocity map (3'.'9±0'.'5 W of 
the nucleus) to the location of the proposed secon d obscured 
mass c oncentration. However, the same authors ijDiaz et al.l 
l2006al . hereafter D06) then performed IFS on a coincident 
but larger region, with better spatial resolution and in the 
near infrared (NIR). They then found that the largest gra- 
dient in the Pa/3 velocity field was in fact some 7" WNW 
of the nucleus and so adopted t his position as a l ikely loca- 
tion for a mass concentration (iDfaz et al]|2006bl . hereafter 
D06b). The most recent work (|Diaz et al.1 12007'') now reports 
three mass concentrations: the visible nucleus and a further 
two hidden mass concentrations at two different locations, 
both apparently more massive than the visible nucleus. De- 
spite the i nconsistenc y between Ha and Pa/3 kinematics (c.f. 
Fig. 5 of iMast et aD and Fig. 3 of D06a), gas kinematics 
in ge neral are well known to exhib it non-gravitational ef- 
fects (jKormendv fc Richstonelll995l ) and the strong velocity 
gradients observed in the gas kinematics could be caused 
by hydrodynamical effects (e.g. outflows, shocks and spiral 
density waves) rather than mass concentrations. 

Nevertheless, the scenario of D06(a,b) is appealing, 
given the theory of dynamical fr iction in a gaseous medium 
as presented by lOstrikerl (|l999') and the numerical imple- 
mentation. Simulatio ns of merging black holes (BHs) in a 
gaseous medium (e.g. lEscala et al.ll2004l : iKim fc Kimll2007l ) 
reproduce density wakes tantalisingly similar to the den- 
sity distributions of the aforementioned nuclear spirals and 
rings; one would be hard pressed to distinguish between Fig. 



5 of lEscala et"aLl (|2004) and Fig. 6 of IMacieiewskil l|2004 ). 
Furthermore, Fig. 2 in lEscala et all l|2004 ) bears a remark- 
able resemblance to the scenario seen at the centre of M83, 
assuming that the high density gaseous wake can initiate 
star formation. As pointed out by D06(a,b), invoking this 
scenario nicely explains the age gradient of HOI (and the 
average age gradient of the HOI data presented in D06b) as 
well as the off-centred old nucleus. 

However, some important queries remain unanswered. 
If the gradient in the gas kinematics is caused by the grav- 
itational potential, one should see the same effect in the 
stellar kinematics which are mostly unaffected by shocks or 
outflows and accurately trace the gravitational potential (al- 
though if the shock and accumulation of gas is significant, 
it may begin to dominate over the smooth potential of the 
disk and bulge). HOI also date some of the youngest clusters 
in the arc to be beyond the present location of the interlop- 
ing mass of D06(a,b). If they are related to the arc of young 
stars, it is not clear how they got to be there: the high den- 
sity w ake c reated by the interloping masses in lEscala et al.l 
('2004') and iKim fc KimI (|2007l '). which we assume leads to 
star formation and an age gradient, is always trailing and 
never preceding the interloping mass. Furthermore, the clus- 
ters nearest the supposed loc ation of the interloper hav e 
been dated at 5 Myrs (HOI and lBresolin fc Kennicuttll2002l '): 
this is larger than one would expect given the age gradient, 
the dynamical timescale and their small distance from the 
interloper (disscussed again in fjG]). 

Perhaps more importantly, there is very little evidence 
of any merger event having occurred recently (the last 50 
Myrs) in M83. T he most promine n t tida l tail of M83 was 
first published bv iMaUn fc HadlevI l|l997h : it is roughly 27' 
in length according to this exposure, subtends around 60° 
in the north by north-west and is estimated to be between 
15' (20kpc) and 20' (26kpc) from the centre of M83. Clearly, 
such an arc is well outside the main galatic disk, defined 
either as the extent of the HII e mission or the Holmb erg ra- 
dius (5^1 and 7'.3 respectively, in lThilker et allbOOSl ). How- 
ev er, this tidal tail is w ithin the HI distribution as seen 
in lPark et all HooJ) and lMifler fc Bregman (2005 '): in fact, 
both these HI stud ies see a no r th-we stern arc of HI emis- 
sion. The figures in lPark et al.l (|200ll ) clearly show this HI 
arc to extend right over from the north-west to the north- 
east and they estimate a radius of 25' from the centre of 
M83. This is ap proximately the same l ocation as the stel- 
lar arc found bv lMalin fc HadlevI (Il997l ) and rescaling and 
overlaying the two images reveals that they are indeed co- 
spatial near the north-west tip of the stellar arc but diverge 
to the north-east with the HI at a larger radius than the stel- 
lar feature. Such features are so far outside the main stellar 
disk of M83 that the disrupted source is unlikely to have 
fallen via dynamical friction to the bottom of the potential, 
but is more likely to contribute to the outer halo. In fact, 
on visualising the aforementioned overlay, one is tempted to 
assign the disturbed outer HI structure seen in iPark et al.l 
(l200ll) to the cap ture of a satellite which we now see as the 
iMalin fc Hadlcv stellar arc. Moving closer to the galactic 
disk, 'Thilkcr ct al. (2005) report GAIA observations of two 
bisymmetric filaments or arms just outside the main HII 
disk (5!l), one extending roughly north from the east of the 
disk and the other extending roughly south from the west 
of the disk, both approximately 15' in length. Such beauty- 
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fully symmetric filaments of star formation are unlikely to 
be induced by a merger event, but rather by the remaining 
HI envelope as the authors propose. The main stellar disk 
of M83 is remarkably regular (as seen in Fig. [l]) with only 
one po tent ial blemish (see j|8l): r adio studies bv lCowan et al.l 
l|l994 ) and lMaddox et al.l (| 20061 ) identify a linear sequence of 
3 sources to the north- east of the central region. However, 
m iMaddox etahl (12006'), the central source of the three is 
identified to be co incident with one of the x-ray sources of 
ISoria fc "w3 (|2003l ) and the authors conclude that this lin- 
ear structure is a background radio galaxy with two lobes of 
emission from jets and a c entral galaxy sourceH 

Finally, the CO map of Sakamoto et al] (|2004l ) is regular 
and symmetric both in its density and its velocity, which is 
unlikely to be the case if a recent merger had disturbed the 
gas sufficiently to ignite star formation in its wake. 



1.6 This Study 

In this paper we report on the results of combining EST 
photometry and ESO VLT NIR long-slit spectroscopy at the 
centre of M83. We address the claims of TOO and D06(a,b) 
regarding the presence of a second obscured mass concen- 
tration by looking for further dynamical signatures in the 
stellar kinematics at the same locations. In particular, to 
help identify a cause for the gradient in the gas kinemat- 
ics seen in D06a, we observe the NIR stellar kinematics at 
the same location. We also build on the study of HOI and 
compare optical and NIR indices of clusters in the circum- 
nuclear arc (some of which were undetected in HOI) with 
models for the stellar populations. Although HOI showed 
that an age gradient existed along the arc from north-west 
to south-east with the youngest clusters in the north-west, 
the reddening vector in this study parallels the tracks in the 
two-colour diagram for ages of 5-10 Myrs, complicating the 
age analysis (Ryd er ct al. 2005. ) , particularly because the 
extinction estimates were derived from the H/3:Ha decre- 
ment which suffers from low H/3 fiux, poor penetration and 
a small wavelength range. NIR indices can probe deeper into 
potentially obscured clusters, giving a more representative 
measure of the ages. Furthermore, clusters can be heavily ex- 
tincted and therefore overlooked in the visible, which is par- 
ticularly relevant for very young clusters (< 5 Myrs) which 
may be dustier and only have a weak intrinsic stellar contin- 
uum (|Rvder et al.1120051 ): this is less of a problem in the NIR, 
but we also present and correct our data using an extinction 
map derived from the Ha:Paa decrement, benefitting from 
more flux, a larger wavelength range and the penetrating 
Paa data. 

The structure of this paper is as follows: 33 describes 
the reduction and homogenisation of data from different in- 
struments; ^describes the analysis of the data; we briefly 
discuss some important details regarding the error analysis 
in 331 SSI presents the results which are subsequently dis- 
cussed in 36] while 33 concludes. 



^ Is is worth noting that there are two other sources roughly 
along the same linear structure in the south-east, sources 32 and 
36, which both have x-ray counterp arts leading to them being 
identified as x-ray binaries (XBRs) in lMaddox et all s study. 



Date 


ID 


PA (°) 


Time 


DIMM 


Data 


D/M/Y 






(mins) 


FWHM 


FWHM 


11/07/00 


[A] 


126.95 


40 (5) 


1'.'02 


0"6 


13/07/00 


[A] 


126.95 


40 (5) 


0'.'52 


0"4 


17/08/00 


[B] 


124.89 


30 (5) 


0'.'74 


0"5 


18/08/00 


[C] 


118.13 


30 (5) 


0'.'61 


0"5 


19/08/00 


[D] 


171.26 


30 (5) 


2'.' 22 


1"2 


21/03/00 


[E] 


51.15 


15 (5) 


0'.'71 


0"4 



Table 1. A summary of the VLT/ISAAC K-band long-slit spec- 
troscopy. The time column gives the total exposure time of all 
frames and, in parenthesis, the exposure of the individual frames. 
The last two columns compare the mean seeing reported by the 
DIMM and the seeing estimated directly from the data f i|2.5ll . 
The ESO program IDs for the March, July and August exposures 
are 64.N-0100(B), 265.B-5723(A) and 65.O-0612(A), respectively 

2 OBSERVATIONS AND DATA REDUCTION 

We present unpublished VLT/ISAAC data and archival 
HST/NICMOS and HST/WFPC2 data. The observations 
and data reduction for each instrument are discussed sep- 
arately below, followed by details of how they were ho- 
mogenised. Fig. [3] illustrates the slit positions relative to 
the HST data. 

2.1 VLT/ISAAC Spectroscopy 

New K band long-slit spectroscopy along 4 different posi- 
tions (A, B, C and D in Fig. [3} was observed in late 2000 
and we combine this with the original major axis data of 
TOO (E in Fig. O to give a total of 5 slit positions in the 
central region of M83 observed with the VLT/ISAAC spec- 
trograph using the 0'.'6 slit (0'.'147 pix~^ plate scale) and an 
SK fllter to give a wavelength range of 2.239vim - 2.362^im at 
a resolution (A/ A A) of 4400. Observations were made using 
the ABBA technique, removing the need for separate sky 
exposures. Table [T] summarises other details regarding the 
exposures. 

A number of kinematic templates were observed with 
the same instrumental setup on previous programmes (ESO 
64.N-0100 and 65.N-0577). However, this template hbrary 
lacked supergiants which are typical of a starburst popu- 
lation. Consequently, we supplemented the library with a 
number of supergiant observations also made with ISAAC 
but using the 0'.'3 slit (ESO 68.B-0530, 67.B-0504). To ho- 
mogenise the library, the latter (of higher spectroscopic res- 
olution) were convolved down to match the resolution of the 
0'.'6 slit templates. The 0'.'3 slit templates also had a shifted 
wavelength range (2.249vim- 2.373vim), so both the 0'.'3 and 
0'.'6 slit templates were cropped to the common wavelength 
range. Our final library consisted of 6 giants (K1-M6) and 
14 supergiants (K5-M5). 

Reduction of the long slit data mainly followed the 
standard technique, as detailed in the VLT/ISAAC user- 
manua0 with use of IRAF0 and ECLIPSEi3 routines. Minor 
differences to the above are summarised below. 

^ http:/ /www. eso.org/instruments/isaac/doc/ 
^ http://iraf.noao.edu/ 

^ http:/ /www. eso.org/projects/aot/eclipse/ 
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Figure 1. A BVR false-colour image of the spiral galaxy M83 (using logarithmic intensity scaling), with the footprint of the ho- 
mogcnised HST/WFPC2 and HST/NICMOS images overlaid in black. The images were observed using the ESO Danish 1.54m telescope 
IILarsen fc Richtlei]|l999h 



The raw data frames were corrected for the odd-even 
effect prior to any other processing using the ECLIPSE Jit- 
ter routine. Persistent bad pixels were detected in the flat 
field images and interpolated over during flat fielding. Subse- 
quently, median filtering was used to detect (but not replace) 
cosmic rays to make a mask for each data frame. After recti- 
fication of the data to a uniform wavelength and spatial grid, 
the frames were aligned and stacked by fitting a Gaussian to 
a star cluster in the flux profile of each frame. Crucially, the 
previous masks were also rectified and shifted in the same 
manner to identify contamination into neighbouring pixels. 



All contaminated pixels were then excluded when combining 
the stacked data frames. Additionally, the data underwent 
an iterative clipping (at 5a) as a precaution to eliminate 
unidentified bad/hot pixels. It is possible to simply remove 
cosmic rays prior to rectification and alignment by interpo- 
lation of neighbouring pixels, but the method employed here 
interpolates from additional realisations of the same pixel, 
which we prefer wherever possible. 

Telluric standard stars (HD119013 and HD118187) were 
also observed each night with the same spectrograph setup 
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Date 


Filter 


Time 


Mean A 


Bandpass 


D/M/Y 






(A) 


(A) 


WFPC2 (Proposal ID: 


8234, PI: Calzetti) 




25/04/2000 


F300W (U) 


700 


3014 


858.3 


02/05/2000 


F487N (H/3) 


1100 


4866 


33.9 


25/04/2000 


F547M (V) 


310 


5488 


637.9 


02/05/2000 


F656N (Ha) 


600 


6564 


28.3 


25/04/2000 


F814W (I) 


237 


8023 


1472.8 


NICMOS (Proposal ID: 7218, PI: Rieke) 


16/05/1998 


F187N (Paa) 


160 


18740 


194.8 


16/05/1998 


F190N (cent.) 


160 


19003 


184.0 


16/05/1998 


F222M (K) 


176 


22181 


1479.4 



Table 2. A summary of the HST observations. Mean wavelengths 
and bandpasses are calculated with bandpar in the STSDAS syn- 
phot package (AVGWV and RECTW, respectively). 



to correct the other observations for at mospheric absorptio n, 
using standard correction techniques (jOrigUa et al.l[l993l ). 

We found that the absolute wavelength calibration of 
the observations varied each night and was not corrected for 
by the ARC exposures. To obtain an accurate calibration 
(the importance of which is described in H2.3|l , we used the 
OH emission line at Avac22460.264A in the unsubtracted 
frames as a reference point and corrected accordingly. We 
further applied a heliocentric velocity correction for each 
night to correct to an absolute velocity scale. 

Slit A was observed over two nights (Table [l]) but the 
seeing varied significantly between them. For all subsequent 
analysis, we chose to use the data with the best seeing (that 
of 13/07/00); the increase in the signal-to-noise ratio (SNR) 
from adding an additional night was marginal and signif- 
icantly degraded the spatial resolution (increasing cluster 
contamination). Problems with the arc lamp observations 
(only the Argon lamp and not the Neon lamp was operating) 
meant the wavelength calibration was also not as accurate 
for the 11/07/00 observations of slit A, further motivating 
our decision to neglect them. The position of slit B, although 
very similar to that of A, covers a slightly difi'erent region 
and is better centred on a few knots, as can be seen in Fig. 
[3l We therefore do not combine the data with Slit A, but 
analyse it separately. 

2.2 HST/NICMOS Images 

The central 20" x 20" of M83 were observed with 
HST/NICMOS in the narrow band fihers F187N (Paa), 
F190N (Paa continuum) and F222M (approximately K- 
band), the details of which are shown in Table [2] Fig. [l] 
illus trates the NICMOS foot print on a false-colour compos- 
ite (|Larsen fc Richtleij[l999l . observed with the ESO Dan- 
ish 1.54m telescope and available on NED). All exposures 
were made with the NICMOS 2 camera with an approximate 
plate scale of 0'.'075 per pixel (see >i2.4p . The data were re- 
duced using the STSCI calibration pipeline. However, for ex- 
posures observed in non-chopping mode (i.e. without back- 
ground exposures), the latter stage of the pipeline (calnicb) 
assumes a sparsely populated field and is configured to es- 
timate and subtract off the background level from the same 



Vrec 
(km s"-') 


Ha Trans. 
Ratio (%) 


[Nil] 
Ha 


H/3 Trans. 
Ratio (%) 


F656N 
F487N 


500 


82.6 


24.6 


96.8 


0.83 


512 


80.6 


25.2 


96.5 


1.03 


550 


72.8 


27.9 


95.8 


1.12 


600 


59.5 


34.0 


94.8 


1.32 



Table 3. The effect of [Nil] contamination and filter transmission 
with recession velocity Vrec for a Gaussian Ha emission line of 
150km s~^ FWHM. The columns represent: the recession velocity 
of the gas Vrec; the Ha transmission at Vrec as a percentage of 
the transmission at Vrec = 0km s~^; the [Nil] contamination as 
a percentage of the Ha emission; the H/3 transmission at Vrec as 
a percentage of the transmission at Vrec = 0km s~^; the ratio of 
the observed (uncorrected) fluxes in the F656N (Ha) filter and 
the F487N (H/3) filters with Vrec (note the effect this may have 
on extinction estimates using these filters). 

exposure. It does this by identifying and masking sources 
(groups of 2 or more pixels 4.5a higher than the image 
median) and using the remaining sigma clipped median as 
the background level. It is not difficult to see that for our 
non-chopped observations of the centre of M83, this process 
will lead to gross over subtraction of the background. Con- 
sequently, the value estimated for the background by the 
pipeline is added back onto the images. They are then flux 
calibrated using the photometric zero-point PHOTFLAM 
calculated by the STSCI pipeline and RECTW bandpasses 
(for narrow-band images) from the STSDAS bandpar rou- 
tine. The background of each mosaiced exposure is estimated 
from a suitable blank region near the edge of the frame which 
is then subtracted. 



2.3 HST/WFPC2 Images 

The central 36" x 36" of M83 was observed with WFPC2 in 
the fihers F814W, F547M, F300W, F487N (H/3) and F656N 
(Ha), details of which are given in Table [2] 

The data were reduced using the STSCI calibration 
pipeline. In addition, the STSDAS routine wfixup was used 
to interpolate over known bad pixels. Flux calibration was 
performed using the PHOTFLAM photometric zero point 
calculated by the STSCI pipeline and RECTW bandpasses 
(for narrow-band images) from the STSDAS bandpar rou- 
tine. 

The transmission curves of the F656N and F814W fil- 
ters were obtained from th e STSDAS calcband rout i ne. Us - 
ing the optical spectra of IStorchi-Bergmann et ahl (|l995h . 
the [Nil] and Ha emission lines were fitted (with Gaus- 
sians) and subtracted; subsequently applying the F656N and 
F814W filter curves gave the ratio of the continuum flux in 
each of the filters. Background subtraction for the F656N 
image was then performed using the F814W image multi- 
plied by this normalisation factor (as done by HOI). 

There are two systematic errors associated with the 
background subtracted F656N observations which required 
correction. Firstly, the F656N fllter is broad enough to in- 
clude [Nil] emission in addition to the Ha emission. Sec- 
ondly, because of the recession velocity of M83, Ha lies near 
the edge of the bandpass, reducing the observed Ha. These 
two systematic errors are also correlated: the recession ve- 
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locity of M83 shifts the Avac6585.27A [Nil] hne out of the 
bandpass, significantly reducing the contribution from this 
line. However, it also shifts the Avac6549.86A [Nil] line into 
the bandpass, increasing the contribution. 

Correcting fo r these two effects is not t rivial . Using the 
optical spectra of IStorchi-Bergmann et all (|l995h . one may 
attempt to correct for both, but the spectra are low res- 
olution (FWHM loA - 500km s"^) causing the lines 
to be significantly broadened. Such instrumental broaden- 
ing hinders our calculations as it is not possible to reli- 
ably estimate the change in transmission for a recession 
velocity of 500km s~^ when the line itself is 500km s~^ 
in width. Therefore we fit the Ha and [Nil] lines in the 
IStorchi-Bergmann et al.l l| 19951 ) spectra with Gaussians to 
estimate their relative fluxes and then reproduce the spec- 
trum with lines in the same fiux ratios but each with a 
FWHM=150km s~^ (in agreement with the ionised gas kine- 
matics of D06). Using this (continuum subtracted) spec- 
trum, we calculate the contribution of the [Nil] lines and 
the transmission of the Ha line in the F656N image for a 
range of recession velocities. Table [3] summarises our find- 
ings. 

Although DOG quote absolute ionised gas velocities 
ranging between 550km s~^ and 650km , the stellar ve- 
locity field is not given, so we have no knowledge of whether 
the ionised gas is actually at a difi'erent velocity to the stars 
or if the absolute wavelength scale is in co nflict with the re- 
cent d etermination of Vrec = 512km s~^bv lKoribalski et al.l 
This latter value is in good agreement with the av- 
erage recession velocity found in the HyperLEDA database 
(507km s~^). We took considerable pains to calibrate our 
ISAAC data to an absolute heliocentric velocity scale (5j2TT|. 
The error weighted mean stellar velocity along all slits was 
found to be 520km s~^. However, the error weighted mean 
molecular Hydrogen velocity w as found to be 5 1 2km s~^. 
We therefore adopt the value of iKoribalski et al.l (|2004l ) for 
the ionised gas and correct accordingly. Quite by chance, the 
[Nil] contamination and bandpass edge effects almost coun- 
teract each other at this velocity, leaving only an overall 
scaling of 0.99 for t he F656N image. 

The spectra of IStorchi-Bergmann et al] l| 19951 ) are av- 
erages over the central region and although we attempt to 
account for the non-uniform gas velocity across the FoV with 
a systematic error (see H4.3p . we might also expect variations 
in the [Nil] :Ha fiux ratio. Thus our correction in this matter 
is only to first order. 

The effect of recession velocity is smaller for the 
WFPC2 F487N (H/3) fiher (Table Ejl and is non-existent 
for the NICMOS F187N (Paa) filter so the F487N image 
received a minor flux correction of 1.036 in accordance with 
Table |3] (neither filter is contaminated by other emission 
lines) . 

2.4 Merging HST/WFPC2 data with 
HST/NICMOS data 

HST narrow band images exist in Ha (WFPC2) and Paa 
(NICMOS). In order to homogenise these two data sets, it 
was necessary to match the pixel sampling, point spread 
function (PSF), position angle (PA) and alignment. Fig. [T] 
illustrates the footprint of the homogenised HST images on 
a false-colour composite. 



A well known feature of NICMOS detector is the rect- 
angular and time variable plate scale. However, the feature 
is well documented and we interpolated the plate scale at the 
time of our observations from the plate scale records avail- 
able on the NICMOS STSCI website. At the time of the Paa 
observations, the plate scale was 0'.'0759788 ± 0'.'0000025 
and 0'.'0752962 ± 0'.'0000025 in x and y, respectively (we as- 
sume that the plate scale of the WFPC2 det ector is constant 
and adopt the value of 0'.' 04554 quoted in iHoltzman et al] 
Il995h . We re-sampled the original 261x61 NICMOS images 
to uniform 261x61 arrays with a 0'.'0752962 pix~^ plate scale 
in X and y. Flux was conserved throughout. 

The WFPC2 images were rotated to the same PA of the 
NICMOS images (using the ORIENTAT keyword) and con- 
volved with a Gaussian PSF of a = 1.355 pixels (the FWHM 
of a diffraction limited PSF for a 2.4m aperture at 1.876n.m 
is 1.503 pixels; if we assume that this Airy disk can be ap- 
proximated by a Gaussian, by adding cts in quadrature we 
reproduce the FWHM of the NICMOS F187N PSF). The ro- 
tated and convolved 800x800 planetary camera images were 
then re-sampled to the same plate scale as the re-sampled 
NICMOS images. 

Finally, we shifted and aligned the rescaled WFPC2 
F814W image to the NICMOS F222M image by maximis- 
ing the cross-correlation function. We applied the same shift 
to the other WFPC2 images and extracted 261x61 images 
corresponding to the NICMOS FoV. The end result was a 
remarkably accurate scaling and aligning of WFPC2 and 
NICMOS data, demonstrated in Fig. [2] The resolution of 
the final images (e.g. Figs. [H [3] Hand [7]) is set by the NIC- 
MOS instrument and we estimate the FWHM to be around 
two pixels (0'.'15). 



2.5 Synthetic Slits: merging space and ground 
based data 

In order to merge the ISAAC data with the NICMOS and 
WFPC2 datasets, we extracted apertures from the HST data 
corresponding to the ISAAC slit positions. As the NICMOS 
and WFPC2 datasets were already homogeneous ( t]2.4p . it 
remained to reproduce the seeing conditions and locate and 
extract synthetic slits. 

To facilitate this procedure, a routine was written to 
minimise the difference between the ISAAC long-slit fiux 
profile (short-ward of the CO bandhead at 2.3/im) and a 
fiux profile produced by overlaying a 0'.'6 slit aperture on 
the NICMOS F222M image, incorporating sub-pixel shift- 
ing, rotation and Gaussian convolution to emulate the seeing 
conditions. The problem was highly non-linear, with many 
local minima which could stall the process. Thus the op- 
timisation required a good initial guess for the parameters 
(made by eye). 

The principal result of this routine was the extraction 
of series of synthetic slits from the HST images, aligned 
to the ISAAC data and matched in spatial resolution. A 
by-product was an accurate estimate for the seeing of the 
ISAAC ground-based data (Table[l| which allows for greater 
accuracy and less contamination when extracting spectra of 
individual star clusters ( ^3.11) . The reported seeing from the 
Differential Image Motion Monitor (DIMM) did not always 
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Figure 2. The homogenised WFPC2 F814W image, shown with logarithmic colour map scaling. Over-laid are contours of the ho- 
mogenised NICMOS F222M image at levels (1.3 X IQ-^^, 1.8 x lO-^^^ 2.5 X IQ-^^^ 7 4 x IQ-iS). The axes are scaled in pc and all fluxes 
are given in ergs s"'^ cm~^. Note the accuracy of the spatial scaling and alignment between the WFPC2 and NICMOS data. The 
footprint of this image is shown in Fig. [T] 

represent the true seeing of our data, which is not surpris- 
ing given that the DIMM looks in a roughly fixed direction 
l|Sarazin fc Roddierlll990l . less than 30deg from the zenith) 
and at visible wavelengths; the ISAAC data was taken at a 
different position on the sky, with a different airmass and in 
the NIR (recall that seeing cx A""'^). 



3.1 VLT/ISAAC data 

The ISAAC K-band long-slit spectra hold both kinematic 
and chemical abundance information. We discuss the ex- 
traction of this information from the ISAAC spectra after 
presenting details of our apertures. 



3 DATA ANALYSIS 

We wish to build on the study of HOI and compare optical 
and NIR indices for clusters in the circumnuclear arc with 
stellar population models and also address the claims of TOO 
and D06 regarding the presence of a second obscured mass 
concentration. The first of these goals requires chemical and 
kinematic analysis of individual clusters. The latter requires 
kinematic analysis along Slit A (or B): Fig. [3] illustrates the 
positions of the putative hidden mass concentrations and 
our ISAAC long-slit data. 

We discuss such analysis below, in two sections: one for 
ground based VLT ISAAC spectra and one for the space- 
based imaging. 



3.1.1 Cluster Apertures 

To study the properties of the individual star forming clus- 
ters at the centre of M83, individual spectra were extracted 
from the long-slit ISAAC spectra. To maximise the SNR and 
minimise contamination from background stars and nearby 
knots, apertures were centred on the peak flux with widths 
equal 2a (!j23}. However, the poorer seeing-limited resolu- 
tion of the ground-based ISAAC spectra sometimes resulted 
in the blending of clusters and the corresponding spectra. 
This is made clear in Table |4] for clusters that are resolved 
into multiple components at HST resolution. Most of the 
clusters seen in the K-band spectra are clearly coincident 
with those detected in HOI. However, some are not; they 
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were too heavily obscured by dust to be detected in the vis- 
ible. Table |4] identifies clusters already identified in HOI. 



to the mean result and then applied a correction using the 
measured stellar velocity dispersion of the cluster. However, 
for the most part, this has negligible effect on W[CO] due 
to the very low dispersions of the star clusters. 



3.1.2 Stellar Kinematics 

The stellar kinematics were extracted both for individual 
clusters and along the who le slit using a Direct Pixel Fitting 
technique ( based on that oflRix fc Whitelll99^ and described 
in detail in lHoughton et aLlbOoi K 

The library of stellar templates used to extract the kine- 
matics contained 20 stars: 6 giants (KlIII to M6HI) and 14 
supergiants (K4I to MSI). To accurately extract kinematics 
of a young population (flux dominated by supergiants), it 
is essential to have a large number and range of supergiant 
templates to minimise template mismatch. This is particu- 
larly important for the CO ro-vibrational absorption after 
2.3ti.m because the depth of the features is strongly depen- 
dent on the spectral type and there is more intrinsic scat- 
ter in the sup ergiant relation compare d to that of giants 
l|Kleinmann fc Hall 1986; Origliaetal] [l993). The stellar 
kinematics were extracted in the range 2.254jxm- 2.355(i.m 
(a limit imposed by the stellar templates observed with the 
0'.'3 slit as described in tj2.1l and by the need for good data 
away from the detector edge). 

When extracting stellar kinematics along the entire slit, 
it was sometimes necessary to bin-up the data; we chose to 
bin the spectra to a minimum SNR of 20 (in the continuum, 
just before the CO bandhead at 2.295jJ.m). 



3.1.3 Gas Kinematics 

The neutral gas kin ematics of the H 2 2-1 S(l) emission 
line (Ava,c22477.17A. lHinkle et al.|[2000h were also extracted 
along the slits. This emission does not coincide with the 
wavelength range used to extract the stellar kinema tics, pre- 
venti ng us from extracting them simultaneously (S arzi et al.l 
I2OO6I ). We instead fitted a Gaussian profile and linear con- 
tinuum to the range 2.243M.m - 2.258^m (-500km s"^ to 
+ 1500km s~^ in velocity; approximately ilOOOkm s~^ from 
the average gas velocity). 



3.1.4 WfCOJ 

We measured the equivalent width of the first CO band 
(^ = 2 ^ 0), W[CO] (as defined bv lOriglia et~aLl Il993l ) 
from the ISAAC K-band spectra. W[CO] was calculated in 
the traditional sense, 



W = 



- 1 



5A, 



(1) 



where fi is the observed (line and continuum) fiux of each 
bin, Ci is the (estimated) continuum fiux of each bin, 5X is the 
width of each bin (in A) and i indexes a suitable wavelength 



range. 



However, lOliva et al.1 (|l995l ) showed that the measured 
W[CO] of a galaxy is affected by the dispersion of the sys- 
tem. We corrected for this using the same technique as in 
iHoughton et al.l l|2006l ): we artificially broadened the stel- 
lar kinematic templates, fitted a second order polynomial 



3.2 HST/NICMOS and HST/WFPC2 data 

3.2.1 Extinction map 

Using the homogenised Ha (WFPC2) and Paa (NICMOS) 
images, we were able to calculate a deep extinction map 
for the centre of M83 as well as for individu al clusters. We 
used a standard Milky Way e xtinction curve (|Cardelli et al.l 
I1989D . recombination ratios o f lOsterbrockl (Il989l ) and a ratio 
of total to selective extinction, Ry, of 3.1 to calculate the 
extinction Ay. Thus we find: 



A^"" =3.1 E(B 



Fh«/Ff 



.13 



(2) 



where Fua and Fpaa are the line fluxes of the Ha and 
Paa emission, respectively. The superscript 'gas' identifies 
that this term has been calculated from the ionised gas. 
Calzetti et al.l (119941 ) found a substantial difference in opti- 
cal depths between the Ha and H/3 Balmer emission lines 
and the stellar continuum underlying these two Balmer lines 
(in starburst galaxies). This was interpreted as a conse- 
quence of the hot ionising stars being associated with more 
dustier regions than the bulk stellar population. Our A'^" , 
determined from the Ha:Paa decrement, will suffer a similar 
effect and we interpret it as the extinction suffered by the 
ionised gas, not the extinction of the (continuum producing) 
stellar population. For brevity, we now drop the superscript 
gas in Ay and E{B — V). 

Figure |4] illustrates the extinction map. A substantial 
fraction of the pixels in the map have no secure upper limit; 
that is to say, within 2a errors, they are unconstrained. A 
similar result evident in Fig. |8] 



3.2.2 Cluster Apertures 

Apertures with the exact same size and position as those for 
the ISAAC data (iJSTTTl} were extracted from the synthetic 
HST slits ( i]2.5|) . to calculate the equivalent widths of the 
Paa and Ha emission for individual clusters. 

As W[Ha] and W[Paa] were derived from homogenised 
narrow band images, we calculate an approximation to the 
true equivalent width, namely 



W: 



AA 



(3) 



where Fx is the observed flux density (line and continuum) 
of the narrow band image, C\ is the flux density of the 
continuum and AA is the bandwidth of the narrow band 
fllter. 

These indices also require correction. As mentioned 
previo usly, when studying starburst galaxies, ICalzetti et al.l 
(|l994l ) found different optical depths for the Ha and H/3 
Balmer emission lines and the continuum underlying these 
two Balmer lines. Consequently, the observed equivalent 
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Figure 3. The 5 slit positions overlaid on the F222M NICMOS image. Positions of the putative hidden mass concentrations are also 
shown as purple (TOO) and red (D06b) triangles. For the latter position, two triangles are plotted: the upright triangle gives the position 
determined from the astrometry while the inverted triangle gives the position estimated directly from the figures of D06b (see i|4.6l l. Flux 
values are given in ergs s~^ cm~^. The position, scaling and orientation is the same as Fig. |2] 



widths of ionised gas emission are affected by this differen- 
tial extinction between gas and stars and require correction 
before comparing to models. 

Using the prescription described in lCalzettil (|l997l ) and 
ICalzettil (|200ll ). the extinction corrections for W[Ha] and 
W[PaQ] were calculated to be 



log 



log 



W(Hq)o 



W(Ha), 



W(PaQ)o 



W(PaQ) 



-0.4E(B-1/) 



= -0.04 E(B - V) 



(4) 



(5) 



where Wo is the observed equivalent width, Wi is the in- 
trinsic width and E(_B — V) is calculated from the HaiPao; 
decrement using a 'standard' extinction curve (333111 ■ The 
errors in the corrected values were nearly always dominated 
by the error in E{B — V). 



3.3 Starburst99 Ages 

SB99 models enable us to age-date individual knots seen 
in the star forming region of M83. We use three indices 
for this purpose: W[CO] direct from the stellar population 
and W[Ha] and W[Paa] from the gaseous nebula emission 
(which essentially measure the same quantity, but differ due 
to the effects of ext inction). W[CO] is defined according to 
lOrigha et all (Il993l ) in both SB99 and our measurements. 
Although W[CO] is affected by the spectral resolution of 
the observations, the low resolution of the SB99 model spec- 
tra in the NIR is not an issue because SB99 calculates 
W[CO] from expan ded versions of the (higher resolution) 
lOriglia et al.l l| 19991 ) models, whic h in turn implemen t the 
theoretical equivalent widths from lOriglia et aU l|l993l ). 

Although many nebula emission indices are calculated 
by default in the SB99 mode ls, W[Paa] is not. Therefore , 
we follow the prescription of iLeitherer fc HeckmanI (Il995l . 
and references therein) and calculate it directly from the 
number of ionising photons (with wavelength < 912A) and 
the model Paa continuum at 1.876^m as a function of time. 
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Figure 4. The extinction Av calculated from the HOfiPaQ; ratio. 
The position, scaling and orientation is the same as Fig. |2] Only 
pixels with 2a detections in either the HQ! or Patt fluxes arc 
shown (pixels not matching this criterion are coloured black). A 
substantial fraction of these pixels have no secure upper bound 
on the extinction value, within the 2a errors; almost none have a 
secure upper limit within 3o". This is a result of the low SNR of 
the NICMOS Paa data. 

3.3.1 SSP models 

Simple instantaneous single stellar population (SSP) burst 
models for WfPaa] and W[CO] (Fig. O dashed lines) were 
calculated with a Salpeter IMF (a — 2.35), a mass range 
of l Mp,< M^. < lOOMp ) and for a metallicity Z = 2Zq 



l|Kobulnickv et al" 19991. appropri ate for M83) using the 



Geneva tracks l|Mevnet et al.lll994 ). No W[CO] is expected 



before 6 Myrs, which is the point at which the massive stars 
of a SSP first evolve to the red super-giant (RSG) phase; 
this property is relatively insensitive to variations in Mup be- 
tween 30 and 100 M0 because of the small nu mber of RSGs 
evolv ing from progenitor masses above 30 Mq l|Origlia et al.l 
Il999t) . These SSP models were found to be incompatible 
with our data for individual clusters: W[Paa] gave consis- 
tently younger ages compared to W[CO], by order of a few 
Myrs. If the RSG phase was occurring earlier than predicted 
by the models, one would expect all W[CO] measurements 
to be > 15 as the time period for which young stars have 
< W[CO] < 15 is very small. We therefore inferred that 
the variation in W[CO] comes from a mixed population, of 
stars that have evolved to the RSG and stars that have not, 
which in turn implies a mix of ages and a finite formation 
timescale for the clusters. 



3.3.2 Mixed population models 

SB99 models only provide two cases of the star formation 
law: an instantaneous burst or constant star formation. We 
therefore investigated the effect of different formation sce- 
narios using the first of these limiting cases. In particular 
we investigated the evolution of a short, finite episode of 
star formation. Similar star formati on laws (albeit for expo- 
nentia l bursts) were investigated in lForster Schreiber et al] 
(|2003l ). but using different models and data. 



All indices for the mixed population were calculated by 
convolving SSP model fiuxes (or flux densities) over time 
with a suitable kernel. Convolution with a top hat kernel 
thus mimicked a finite episode of constant star formation. 
W[PaQ] was calculated by convolving the number of ionis- 
ing photons of a SSP and the related model continuum flux 
density at 1.876/im with the desired kernel. We then cal- 
culate W[PaQf] as before using these new mixed population 
values. When calculating W[CO], we first converted it into a 
flux using the model continuum at 2.29/^m. This wavelength- 
integrated flux and the same SSP continuum at 2.29/im were 
then convolved with the desired kernel and used to recalcu- 
late W[CO] for the mixed population. 

Fig. [5] shows Wpaa] and W[CO] for an episode of con- 
tirmous star formation lasting 6 Myrs (solid lines) . The effect 
of a finite duration star formation on W[PaQ] is clearly neg- 
ligible, but the onset of the W[CO] index is brought forward 
in time by a few Myrs. As expected, after around 10 Myrs, 
the SSP and mixed population models are indistinguishable 
as changes in the indices become slower and less pronounced. 
The mixing timescale of 6 Myrs is somewhat arbitrary and 
was chosen by eye to to best fit the data. Formal errors are 
difficult to estimate, although the data appears to be incon- 
sistent with timescales less than 5 Myrs and larger than 7 
Myrs, giving a 1 Myr margin. 

Strictly speaking, we cannot perfectly reproduce a fi- 
nite duration of constant star formation from instantaneous 
burst models; rather, we are in fact producing a mixed pop- 
ulation via many sequential instantaneous bursts over a fi- 
nite duration and assuming that this closely approximates 
a short episode of constant star formation. The time res- 
olution for the mixing (separation between the sequential 
SSP models) is 0.1 Myrs. We see no significant change with 
smaller time resolutions and so presume that we closely ap- 
proximate the ideal case of continuous star formation over 
a finite time. 

Using the same techniques, we investigated the effect of 
the functional form for the star formation during the finite 
episode and found reasonable fits for a finite burst of expo- 
nentially increasing star formation (e-folding time ~ 3Myrs) 
and a finite burst of exponentially decaying star formation 
(e-folding time ~ 2Myrs). Although the finite episode of con- 
stant star formation always provided an overall better fit, 
our data was not able to reliably distinguish between such 
alternatives, so we chose to adopt the simplest scenario (i.e. 
a short burst of constant star formation). We also tested 
the effect of contaminating an SSP model with an old (Gyr) 
population and such a scenario could not reproduce the ob- 
served data trends: the W[CO] estimates were too high for 
the younger clusters (i.e. those with high W[PaQ:]). Finally, 
we tested multiple SSP bursts separated by the dynami- 
cal timescale and half the dynami cal timescale, motivated 
by the work of Allard et all l|2006l '): [Sar zi et al. (2007) and 
[Falcon-Barroso et al.l ( 2007l l and our attempts mostly suf- 
fered the previous problem: too much W[CO] at younger 
ages. However, although many episodic bursts (i.e. more 
than 2) were unsuccessful at reproducing the observed data 
trends, models with only two burst were more successful: 
fine tuning the mass fraction of the latter burst to be 3 - 4 
times the initial burst mass, setting the delay between the 
bursts to 5 Myrs and diluting the W[CO] with a fraction of 
the total luminosity to mimic hot dust emission (somewhat 
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simp listic compared to the models of lForster Schreiber et al] 
I2OO3) and replacing the SSP bursts with finite (Myr) bursts 
went some way to reproducing the data but did not com- 
pletely resolve the discrepancy. Thus, despite considerable 
fine tuning of various parameters, multiple bursts did not 
match the data as well as a constant burst lasting 6 Myrs. 

SB99 was rec ently updated to off er models based on ei- 
ther the Geneva (iMevnet et al. II1994') or the Padova tracks 
(jFaeotto et al ] 1 1994 ). Ahhough we illustrate models using 
both the Geneva and Padova tracks in Fig. O we use Geneva 
tracks to estimate the ages of the clusters in the final anal- 
ysis (the same as HOI); the Padova tracks with a similar 
metallicity (2.5Z0) still match the data and provide a simi- 
lar age gradient, but with a tendency for younger ages. This 
is explained by the slightly earlier onset of W[CO] in Fig. [5] 
when using the Padova tracks. 

It is not trivial to read the cluster ages from an age- 
metallicity grid such as Fig. |6] if one allows variation in 
the metallicity, there is a larger range of acceptable ages 
given the age-metallicity degeneracy. As the random errors 
in W[Paa] are usually much larger than those of W[CO], 
one could use a fixed metallicity {2Zq) and take the nearest 
age along the horizontal but such an approach is unfeasible 
when the error for W[PaQ:] is small compared to the hori- 
zontal distance to the nearest model track (as is the case for 
cluster 8). We derived ages from models of a fixed metal- 
licity {2Zq) using the shortest line from the observation to 
the 2Z0 model tracks along W[Paa] and W[CO], effectively 
ignoring the larger uncertainty in W[Paa] and any variation 
in metallicitj(f|. We show the Zq models (and how they link 
with the 2Zq models) in Fig. [6] to justify our metallicity 
assumption: most of the clusters scatter around the 2Zq 
tracks. 



4 ERRORS 

The estimation of errors in our datasets deserves special 
mention as they sometimes required care and are not always 
calculated using the standard first-order approach, for good 
reason. 



4.1 Homogenised HST data 

By homogenising the HST datasets, we introduced correla- 
tions between the pixels and between the errors. Whenever 
we added pixels, we added the variances. This would usu- 
ally cause our errors to be severely underestimated given 
that the variances become correlated and would not add 
in such a manner. However, careful scrutiny of our method 
using Monte-Carlo simulations of noise frames actually in- 
dicates that our final random errors are over-estimated by 

^ Perhaps a better method would be to weight the difference in 
W[Paa] and W[CO] between the models and observation with 
the observational errors and then find the closest model point 
using these scaled distances, but in our case, the rewards of the 
increased accuracy would be minimal and somewhat futile con- 
sidering the systematic errors discussed in ^4.51 
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Figure 5. W[Paa] and W[CO] predictions from SB99 with a = 
2.35 and Mup = IOOMq. Dotted lines represent the instantaneous 
burst models. Solid lines represent mixed models for a finite burst 
of constant star formation over 6 Myrs. The time scale on the plot 
is relative to the mid-point (mass we ighted mean) of the burst. 
Models for both G eneva CG lMevnet~ et al. 1994) and Padova (P 
iFagotto et al. 1994) tracks arc shown, thanks to the recent update 
to SB99 i fVazquez fc Leitherer.i2005() . 



a minimum factor of ~ 1.7. The cause is the seeing convo- 
lution stage: we simply convolve our variance arrays in the 
same manner as the data, but this only results in obtaining 
a seeing weighted average variance for each input pixel; in 
reality, the true error is reduced much more than this. Thus, 
the greater the seeing, the larger our over-estimation of the 
errors. We estimate factors of (1.7, 1.9, 1.9, 2.9) for slits (A, 
B, C, D), respectively but we do not apply these corrections. 

One could question why we degraded the Ha and Paa 
data to match the spatial resolution of the ground-based 
data: we could obtain better spatially resolved data for these 
indices. However, then one is faced with the dilemma of com- 
bining non- homogeneous data and which to trust. The final 
spatial resolution of our data is thus limited by the seeing, 
but it is worth recalling the resolution of the study by HOI, 
which is only a factor of two better. Although the diffraction 
limited FWHM of HST at 656nm is 0'.'056, the pixel sam- 
phng of WFPC2 is 0'.' 04554. However, HOI further convolved 
the data with a Gaussian of a = 2 pixels (FWHM = 0'.'21) to 
create sufficiently uniform cluster profiles that could be well 
fitted by a single PSF model. This brought the final spatial 
resolution of the HST data to 0'.'22. The resolution of our 
best data is 0'.'4 — only a factor of two larger. Furthermore, 
HOI calculated W[Ha] within 0'.'56 and l'.'34 diameter aper- 
tures (12pc and 29pc at 4.5Mpc) and found little difference 
between them for most clusters. Our apertures are generally 
within this range: they are rectangular with heights equal 
to the slit width (O'.'e) and widths equal to (O'.'S, 0'.'4, 0'.'4, 
I'.'l) for slits (A, B, C, D), as determined by the seeing. We 
choose to use the higher resolution of the space-based data 
only to flag potentially contaminated (unresolved) data in 
Table |4l Thus we measure luminosity weighted ages for any 
contaminated clusters which may be particularly true for 
clusters 5, 8 and 12 in Table H 
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4.2 W[CO] 

If we calculate the random error in W[CO] from error in 
the flux values alone, we find that it is negligible compared 
to the error introduced by uncertainly in the velocity and 
dispersion. Therefore, we use zblcr random errors in the ve- 
locity and velocity dispersion measurements of each cluster 
to calculate the systematic error in W[CO]. We then use 
whichever is the larger. As a result, the error in W[CO] is 
often asymmetric. 

4.3 Ha 

As discussed in ^2.31 the Ha flux is subject to an uncertainty 
in the gas velocity due to the systematic velocity of M83 
placing the emission on the edge of the F656N bandpass. We 
therefore appoint minimum errors of 10% to all Ha fluxes 
(often dominating the error in W[Ha] for the clusters) and 
in calculating the extinction estimate of Fig. [8] This 10% 
systematic error corrects for velocity range of 430km s~^ 
to 560km s~^. Because this error is systematic, it does not 
change with binning or aperture size and is applied as a last 
step before analysis. 

4.4 Av, E(B-y) 

Extinction estimates using a HH decrement are often tricky 
due to the uncertain penetration depth of the emission lines 
in question: estimating the extinction in dusty regions using 
the H/3:Ha decrement invariably leads to a lower estimate 
than comparing the Ha:Paa decrement as the effective pen- 
etration is much greater with the Paa line. We assume that 
our Ha: Paa decrement is not limited in this way for the 
range considered here (0 <Av< 8). However, estimating the 
random error is still non-trivial. 

A standard first-order approach gives that the error in 
Eq. [2]is CTA^ = — 1.6(Fpaa/FHa). However, this is only valid 
for small errors in Fpaa and Fhq: the errors in Paa data 
are large (sometimes S> 10%) and the approximation breaks 
down. We therefore calculate the error explicitly by substi- 
tution. Due to the non-linear nature of Eq. [21 there is con- 
siderable difference between the la and 2(t confidence limits. 
We therefore opt to quote 2(t values because Icr values are 
significantly smaller and lead to a false impression of accu- 
racy. The errors in the corrected W[Ha] and W[Paa] of Fig. 
[S^c) and (d), dominated by the error in Av were calculated 
using this 2a error. 

4.5 Cluster Ages 

We wish to quote an error for the ages of the clusters, but 
this is not trivial. Here, we list the reasons for the difficulty 
and explain how the age errors in Table |4] are determined. 

The cluster ages are calculated from the NIR indices 
W[PaQ] and W[CO]. The NIR indices (particularly W[Paa]) 
suffer less from extinction and are therefore more accurate in 
terms of random error and also any systematic error associ- 
ated with the extinction correction (e.g. the extinction law, 
the unknown penetration depth of our Ha and Paa data) 
because the correction is smaller. W[CO] is corrected for ve- 
locity dispersion effects but this correction is very small and 
does not affect the derived ages. However, young clusters 



are expected to be dusty and if heated by UV radiation, 
this dust re-emits, mainly in the mid-IR but also in the NIR 
which could diminish the W[CO] index by artificially in- 
creasing the continuum level and thus affect our derived ages 
and the mixing timescale of i]3.3.2l however, this effect is 
usually very weak and we do not att empt to correct our data 
for th is effect, unlike the study bv iForster Schreiber et aU 

(Hooi). 

The choice of evolutionary tracks also plays a role: we 
quote ages using the Geneva tracks, but the Padova tracks 
also fit the data albeit with systematically lower ages of 
around 0.5 Myrs for the youngest clusters (4.5 Myrs old) 
and 1.5 Myrs for the oldest clusters in the arc (8 Myrs old). 
This difference leads us to quote a systematic uncertainty in 
our age estimates of order zblMyr. 

We also vary the mixing timescale of the models to 
match the cluster data ( H3.3.2p . The mixing timescale was 
chosen by eye to fit the observations, which were inconsis- 
tent with mixing timescales less than 5 Myrs or larger than 
7 Myrs. Error in this parameter will manifest itself as a 
systematic error in the ages, particularly for the youngest 
clusters as this is where the change is greatest. We also note 
that without applying the Calzetti et al. extinction law, the 
data of Fig. [Gja) is best fit by a burst duration of 5.5 Myrs. 
Thus, given that the mixing timescale is estimated to be 
accurate to a Myr, we estimate the youngest clusters could 
also suffer another systematic error of order ±lMyrs. 

We are also somewhat limited by our ground-based spa- 
tial resolution. Because we homogenise all our data, this af- 
fects W[Ha], W[PaQ] and W[CO] and the derived ages. If 
two clusters are unresolved, we expect to measure a lumi- 
nosity weighted W[Ha], W[Paa] and W[CO]. But whereas 
younger clusters (< 5 Myrs) will have higher Ha and Paa 
emission (after correcting for extinction), the K-band lumi- 
nosity and the W[CO] is dominated by older (10 Myr) clus- 
ters. Hence if the two unresolved clusters are different in age, 
we obtain conflicting ages from the different indices. How- 
ever, one expects neighbouring clusters to have similar ages 
(given the results of HOI), so this effect is probably small. 
Nevertheless, it is possible to identify contamination from 
another cluster or from the bulge stars (for highly extinct, 
dim clusters). The luminosity weighted stellar velocity dis- 
persion of low mass star clusters is very low (< 20km s~^), 
so if we measure them to be high (> 30km s~^), the K- 
band light is likely not dominated by a single cluster. Con- 
sequently, when a cluster has a large velocity dispersion 
(> 30km s~^), it is noted in Table |4l Furthermore, when 
extracting the kinematics of each cluster, we fit a linear 
combination of stellar templates. As discussed in i|2.1l we 
have a large number of giant and supergiant templates to 
fit the spectra and the best-fit template can be used to in- 
dicate the different stellar populations in the spectrum, and 
the relative fractions. Therefore, in Table |4] we give the frac- 
tion of giant and supergiant stars used to create the best-fit 
kinematic template for each cluster. The giant/supergiant 
fraction clearly correlates with the velocity dispersions and 
we can identify which cluster measurements are contami- 
nated, and to what degree. What effect the fraction of con- 
tamination has on the derived ages is not always obvious, 
but is likely to lead to bias towards the age of the contami- 
nant; thus we quote the ages (from HOI) of any contaminat- 
ing clusters resolved in HST images in Table [4] so that the 
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reader can compare the measured age of the cluster with the 
measured age of the contaminant (s). 

Finally, we have to match the observed cluster data with 
the models. Given the age-metallicity degeneracy, this is not 
trivial and we discuss our method in i)3.3.2l The proximity 
of our observations to the nearest point on the 2Zq model 
track is an indication of how well the model fits the data, but 
it is also an indication of the stochastic nature of star for- 
mation: given a continuous IMF and only a relatively small 
quantity of gas, the stars formed will not be uniformly dis- 
tributed over the entire mass range but will be scattered at 
random positions weighted by the IMF (especially true for 
the red supergiants). Thus we expect an intrinsic scatter in 
the number of red supergiants per cluster given a known 
mass and age, which leads to an intrinsic scatter in the ob- 
servable quantities WfPaa] and W[CO] — they will never 
agree perfectly with the models. Given that HOI estimate 
most of the clusters to have masses around 10''-10'^Mq, we 
must consider these stochastic processes. For the majority of 
cases, the agreement between data and models is very good 
which leads us to estimate a random error of 0.5 Myrs for 
most of the clusters. However, there are some clusters which 
are more difficult, such as 8 and 10; for these clusters we 
do not know if error (systematic or random) or nature has 
resulted in the data falling far from the models. Somewhat 
arbitrarily, we assign random errors of zblMyr to these clus- 
ters. For clusters that do not match the model predictions 
for both W[Paa] and W[CO] (clusters 14 and 15) we choose 
to trust the W[CO] index (see ^5.1|l but assign asymmetric 
errors to cover the range of all possibilities. For clusters that 
are do not match the models and are assumed to be very old 
(> lOMyrs) we do not quote an error, only a lower limit. 

4.6 Astrometry 

Fig. [3] shows the positions of the putative hidden mass con- 
centrations relative to the NICMOS data and the ISAAC 
slits. 

The slits were located on the NICMOS image as de- 
scribed in H2.5I and were not positioned using the astrom- 
etry provided with the ISAAC and NICMOS data. Formal 
uncertainties are difficult to calculate although by trial and 
error, we estimate the error in x and y for Slit A to be at 
around O'.'l (comparable to a NICMOS pixel) and the error 
in the position angle to be at most 1 deg. However, these val- 
ues vary for each slit: with poorer seeing, the ISAAC data is 
smoother and the errors in the slit positions increase slightly. 

The pipeline reduced WFPC2 and NICMOS images 
come with world coordinate systems. However, the two sys- 
tems were not comparable: the cluster positions of HOI were 
correctly centred on the WFPC2 images (as to be expected, 
for they were derived from the same data) but they were off- 
set when overlaid on the raw NICMOS images (0'.'54 to the 
SE as determined from the position of the nucleus in the 
raw F814W image and the raw F222M image). We chose 
the path of least resistance and adopted the WFPC2 as- 
trometry. To calculate a world coordinate system for our 
homogenised data, we used our precise knowledge of the 
pixel size of the homogenised HST data (square pixels with 
sides 0'.'0752963 in length) and varied the absolute refer- 
ence to best align the positions of the clusters from HOI on 
the data. We assumed no error in the position angle of the 



WFPC2 and NICMOS data as defined by the ORIENTAT 
header keyword. Comparing our astrometry, we find the nu- 
cleus at (13h37m00.91s, -29d51m55.7s) whereas D06b find 
it at (13h37m00.95s, -29d51m55.5s): an offset of 0'.'49 to the 
NE. 

With no agreement between the WFPC2, NICMOS or 
D06b astrometry, we retain that of the WFPC2 images. 
However, when plotting the position of the hidden mass con- 
centration given in D06b we apply a shift of 0'.'49 to the SW - 
the same shift which would align the position of the nucleus 
given by D06b to the position of the nucleus in our world 
coordinate system (WCS). Quite by chance, more than by 
design, this places the position of the putative mass con- 
centration securely in the middle of Slit A. However, D06b 
state a 2a error of 0'.'7 in their coordinates for this position; 
our slit width is 0'.'6. Thus, it is reasonable to argue that the 
absolute centre of the proposed mass concentration could 
still lie outside our slit position; however it is unreasonable 
to argue that we wouldn't detect its gravitational influence 
in the stellar kinematics given that the gas kinematics ap- 
pear kinematically disturbed over many arcseconds in data 
of D06. For completeness, in Fig. [3] we also plot the posi- 
tion of the DOB interloper using their images as a reference, 
rather than their coordinates; this second location still over- 
laps that of Slit A but is approximately 0'.'38 (8pc) west of 
the formal position in our WCS. 

The position of the mass concentration proposed by TOO 
was considerably easier to locate: by using the position of 
the second dispersion peak in the original data and then 
locating the position of this slit on the homogenised HST 
data using the same techniques as for the other slits ( H2.5p . 
we were able to locate the second dispersion peak on the 
homogenised HST data without resorting to any WCS. By 
design. Slit A overlays this same position. We assumed that 
the second dispersion peak was centred along the slit width, 
which could be open to speculation. However, as Slit A is 
roughly perpendicular to the original slit of TOO (Slit E in 
Fig. [3)l , we would adequately cover any positional error per- 
pendicular to the original slit. 



5 RESULTS 

We present our findings with regard to the cluster ages and 
the putative hidden mass concentrations, in two appropriate 
sections. 



5.1 Cluster Ages 

Fig. Oa) shows the mixed model predictions for the evolu- 
tion of W[CO] and Wpaa] for solar (Zq) and twice-solar 
(2Zq) metallicity with data points for individual clusters 
over plotted. The same models and clusters are plotted in 
Fig.Etb) for W[CO] and W[Ha]. As W[Ha] and W[Paa] are 
essentially measuring the same quantity (the number of ion- 
ising photons), the main difference between the two figures 
is the effect of extinction. 

Green points indicate less extinct clusters (Av< 2), that 
roughly match the 2Zq mixed population models for W[Paa] 
W[Ha] and W[CO], whereas the orange points indicate more 
extinct clusters (Av> 2), that tend only to match the model 
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predictions for W[Paa] and W[CO] and lie below the pre- 
dictions for W[Ha]. At first, W[Ha] and W[Paa] therefore 
appear to be uncorrelated in Fig. [6l[a,b) but after applying 
Calzetti's differential extinction, the correlation is obvious 
in Fig. [6l[c,d). The duration of the star formation episode 
(6 Myrs) was therefore chosen to fit both the (corrected) 
W[Hq]-W[CO] and W[PaQ]-W[CO] data in Figs. Hc.d). 

Red points indicate clusters that do not match model 
predictions in the (corrected) W[Ha]-W[CO] and the 
W[Paa]-W[CO] plots. Although only models with ages < 10 
Myrs are given in Fig. (6] these points are not consistent 
with older populations either, due to the larger values of 
W[Ha] and W[Paa]. However, two of these clusters (16 and 
17 in Table |4)| have very low W[Ha] and W[Paa] and it is 
therefore important to consider non-photoionised gas emis- 
sion (caused primarily by shocks and their precursors, from 
supernovae and massive stellar winds but also turbulent 
mi xing layers and chan ges in gas temperature as explained 
by ICalzetti et al.l |2004 ) . Although these two clusters are 
not in the non-phot oionised regions of M83 highlighted by 
ICalzetti et al.1 (|2004l ). one cannot discount the possibility of 
non-photoionised emission. As they are both quite distinct 
from the main star forming arc, we conclude that the Ha and 
Paa emission is non-photoionised and thus W[CO] reflects 
a much older population, greater than 10 Myrs. The other 
two red points (14, and 15 in Table |4]) are somewhat puz- 
zling (HOI also found the photometry of these clusters to be 
anomalous): either the nebula emission equivalent widths or 
the W[CO], or both, are low compared to the models. Com- 
paring the ages predicted by each index, we noticed that the 
ages predicted by W[CO] were in good agreement with the 
ages of nearby clusters, but the ages predicted by W[Ha] 
and W[Paa] were not. Therefore we tentatively assign them 
ages based on W[CO] measurements alone (although cluster 
15 is the only cluster entirely consistent with the instanta- 
neous SSP models, it makes little difference to the derived 
age). 

It should also be noted that clusters 5 and 12 are in 
a particularly densely populated region and their ages are 
likely more representative of the local average (of those clus- 
ters listed as contaminants in Table |3} because of significant 
contamination. Similarly, cluster 9 is probably heavily con- 
taminated by cluster 13, although we do register significant 
differences in W[CO] for the two. 

Fig. [7] illustrates the positions of the clusters with their 
ages (rounded to the nearest Myr) from Fig.[6jc), while Ta- 
bleUquotes them to the nearest 0.5 Myrs. Although the ages 
from Fig. [HIc) and Fig. IHtd) agree remarkably well for each 
cluster, we use the age derived from the NIR data (W[Paa]) 
because it is less effected by extinction and is therefore more 
accurate (as discussed in t)4.5p . 

5.2 Hidden Mass concentrations 

Fig. |3] illustrates the positions of the putative hidden mass 
concentrations of D06 and TOO and our long-slit data. Both 
slits A and B show identical kinematic features at the lo- 
cation of D06b, therefore for brevity we show only the the 
kinematics along Slit A in Fig. |8] 

A sharp gradient in the H2 gas velocity Vgas is found in 
the same vicinity that D06 report a gradient in the Pa/3 gas 
velocity (140pc along slit A) and both gradients are in the 



same direction. The gradient in H2 is sustained (i.e. Vgas is 
roughly constant for 80pc to the north). There appears to 
be a gradient in the opposite direction at 60pc along the 
slit, although the line flux is low, the kinematics of Slit B 
disagree and there are luminous star clusters at this location 
which might bias the observations. 

The stellar velocity V* and stellar dispersion a* profiles 
both show synchronous drops in the vicinity of the D06 mass 
concentration (160pc along slit A), but these changes are not 
sustained and V* and a* both recover, although the recovery 
is significantly more gradual on the north side. 

All these features occur approximately at the edge of the 
north-east dust lane, clearly seen by plotting the extinction 
(Av) along the slit in Fig. |8] 

In Fig. [9l we show a position-velocity diagram for the 
stellar kinematics. The best-fit Gaussian line-of-sight veloc- 
ity profiles (VPs) are shown at each position along the slit 
and each is scaled to the same maximum intensity for illus- 
tration purposes (i.e. the VPs are not normalised). Where it 
was necessary to bin the spectra to extract kinematics, the 
VPs in Fig. [9] are stretched accordingly. It appears from this 
figure that the changes in the stellar velocity and dispersion 
are synchronised with the sudden increase in the extinction 
via the dust lane. However, Fig. [9] suggests that the VPs ex- 
tracted at the edge of the dust lane are principally affected 
at the high-velocity wing; the low-velocity side of the VP 
appears unaffected by the presence of the dust. In fact, if 
the kinematic fits are extended to include /13 and h^, we 
find very asymmetric VPs at the location of the dispersion 
drop with /13 ~ 0.25. 

Further along the slit, we see a very small rise in the 
dispersion at the TOO location and the dispersion here is in 
agreement with that of the second dispersion peak in TOO 
of ~ 70km s~^: there is no clear dispersion peak along this 
slit position. 

Towards the south end of slit A (at around 340pc on Fig. 
[8| we see a further rise in V* and a* accompanied by another 
strong gradient in Vgas. The extinction shows a marginal 
increase in roughly the same vicinity (320pc along the slit) 
but thereafter shows no obvious increase, although our errors 
for the extinction are unbound because of the very small Ha 
and Paa fluxes received here. 



6 DISCUSSION 

Fig. Ogives the clearest evidence yet of an age gradient along 
the circumnuclear arc, with the youngest clusters nearest 
(and coincident with) the north-east dust lane. The ages are 
generally in agreement with those found by HOI, but with 
considerably less scatter. As done in D06b, if the HOI data is 
averaged over angular sections, then the age gradient is very 
similar to the one found here, although this work extends 
to the highly extincted clusters in the dust lane and the 
absolute ages are around a Myr older here. 

We also see from the false-colour image of Fig. [7] that 
the star formation does not stop where the north-east dust 
lane intersects the arc, but continues through it and along 
the northern edge: there appear to be many star clusters at 
the young end of the arc, embedded in the dust lane. HOI 
considered the age gradient along the arc to be an indication 
of star formation propagating from the south-east end of 
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Figure 6. The SB99 models for the mixed population of Fig[5]with data points for individual star clusters. Unlike Fig.O two metallicities 
are shown: the top (magenta) curve shows models for Z=2Zq while the lower (magenta) curve shows models for Z=Z0; the vertical blue 
lines link points of constant age and are spaced in 0.5 Myr intervals. Absolute ages are shown in Myrs on the upper (2Zq) curve. The 
dashed (magenta) line gives the model predictions of the unmixed SSP (2Z0). The data points represent measurements for individual 
clusters: green points indicate low extinction clusters consistent with ages < 10 Myrs; yellow points indicate high extinction clusters 
consistent with ages < 10 Myrs; red points indicate clusters that are inconsistent with ages < 10 Myrs. The upper two panels, (a) and 
(b), show the extincted data; the lower two panels show the same plots with W[Ha] and W[Pa«] corrected for differential reddening 
using the Calzetti prescription ( N3.2.2I I. Panel (b) also quotes Avfor individual clusters. Horizontal error bars in the lower (extinction 
corrected) panels, (c) and (d), are often dominated by the uncertainty in the extinction: 2cr errors are shown for all clusters except for 
two red data points: the extinction, and therefore the corrected W[Paa] and W[Hc(], are unconstrained within the 2a errors, due to low 
Ha and Paa flux. 



the arc to the north-west, terminating roughly at the dust 
lane. D06 alternatively proposed that this propagation of 
star formation was a consequence of the dynamical friction 
trailing an interloping hidden mass concentration, such as 
a satellite galaxy (although see i lG-lfl . However, we propose 
that the arc of young star clusters is a result of gas inflow 
along the north-east dust lane creating a nuclear spiral of 
gas and an associated ring of star formation, because we find 
no evidence to support the presence of an interloper at the 
position of D06(a,b) or the position of TOO. 



6.1 No hidden mass 

Previous claims of hidden mass concentrations are incon- 
sistent with our observations for many reasons which we 
outline below, in addition to the concerns raised in i ^l.51 We 
note that the strong gradient in the molecular gas velocity 
is 20pc north-west (left along slit A) of the formal position 
quoted in D06b for the hidden mass concentration. Whether 
this is a systematic error between coordinate systems or a 
genuine offset is unclear, although, as discussed in 8 ^4.61 close 
inspection of the figures in D06b provides a slightly different 
location, 0'.'4 (8pc) west of the formal position, going some 
way to rectifying the difference. 



A large mass concentration of the magnitude discussed 
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Figure 7. A false-colour imago of the centre of M83 made with homogenised NICMOS and WFPC2 data. The position, scaling and 
orientation is the same as Fig.|2] NICMOS F222M is red, WFPC2 F814W is green and WFPC2 F300W is blue. A non-linear grey-scale 
has been employed to enhance low luminosity (heavily extinct) clusters. Note that the circumnuclear arc can be seen to extend into and 
along the dust lane. Cluster ages (derived from Fig. |6] and rounded to the nearest Myr) are over plotted: colours of the ages aid the eye 
in distinguishing the age gradient along the arc; green represents 5 Myrs, yellow represents 6 Myrs, red represents 8 Myrs while white 
represent uncertain ages. The resolution of this image is set by the NICMOS instrument and we estimate a FWHM of O'.'IS. 



in D06(a,b) would produce an observable rise in the velocity 
dispersion — as seen for the nucleus of similar mass (peak 
value a ~75km s~^). We see no rise in the velocity disper- 
sion at the location of D06 or at the location of the gradient 
in the molecular gas velocity. The highest stellar velocity 
dispersion seen at the position of the gradient in molecular 
gas velocity is 33 ± 2km s~^and just 15pc right along slit A 
is a minimum in the stellar velocity dispersion, measured to 
be 15±3kms-^ This minimum in stellar velocity dispersion 
is probably caused by the young (6 Myrs) star cluster 5 in 
Table |3] but as we shall see, it is highly unlikely that this 
cluster is dominating the light from an extincted mass con- 
centration. The minimum in stellar velocity dispersion may 
also be an effect of dust extinction: iBaes et"al] l|2003l ) have 
shown that dust obscuration can lead to asymmetric veloc- 
ity profiles and apparent drops in stellar velocity dispersion. 



although why we see the drop at the location of the shock 
and not at the peak of the extinction (some 50pc apart) is 
far from clear in this scenario. 

At the location of the putative hidden mass concentra- 
tion, the extinction (Ay) is 4±1 (2 sigma errors). This es- 
timate is derived from the HII emission — the same tracer 
used to obtain the rotation gradient in D06(a,b). In the K- 
band, we therefore expect the extinction Ak to be 0.4 mag- 
nitudes, so we will receive around two thirds of the total flux 
emitted from any stars co-spatial with the ionised emission 
(assuming foreground extinction). Yet we do not see any in- 
crease in surface brightness to indicate an increase in stellar 
density, other than cluster 5, which we measure to have a 
low stellar velocity dispersion of ~ 15km s~^ and an age of 
6 Myrs. Thus we conclude that there is no large ( > IO^Mq) 
hidden mass there. 
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Figure 8. Upper panel: the flux profile of slit A. Positions of the 
putative hidden mass concentrations are also shown as purple 
(TOO) and red (D06b) triangles. Middle panel: The W[CO] mea- 
sured in each bin along slit A: red lines indicate maximum sys- 
tematic errors from velocity and dispersion uncertainty; the error 
bars on the points correspond to random error. Lower panel: The 
Kinematics along slit A: the black, red and green lines illustrate 
the stellar velocity (relative to galaxy's mean stellar velocity), 
stellar velocity dispersion and the H2 gas velocity (also relative 
to galaxy's mean stellar velocity), respectively; the grey line shows 
the extinction. Ay calculated from the Ha and PaQ? emission and 
the lighter grey shading illustrates the 2a error. Both panels: The 
PA of the slit is 126.95°, with the northern side on the left. 
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Figure 9. A position velocity diagram of the stellar kinematics 
from Fig. [8] Where it was necessary to bin the data to extract a 
VP, they have been stretched to represent the space over which 
the data was binned. The extinction profile Av of Fig. |8] is over 
plotted in white. Positions of the putative hidden mass concen- 
trations are also shown as purple (TOO) and red (D06b) triangles. 



We, like HOI, observe the youngest star clusters to pre- 
ceded the trajectory of the putative mass concentration. 
This is contrary to Ostriker's theory of dynamical friction 
in a gaseous medium and the subsequent hydrodynamical 
simulations as discussed in i]1.5l These young clusters, as 
an extension of the arc, cannot be explained by invoking an 
interloping mass. 

Shocks are well known to give rise to radio emission, as 
discussed in i]1.2l and the MIR is well known to be correlated 
with the radio; thus the peaks observed in both wavebands 
near the position of the putative interloper are easily ex- 
plained by a shock. 





200 300 
Dist. (pc) 

Figure 10. As for Fig.|8]but for slit E - the original data of TOO. 
The position of the putative hidden mass concentration from TOO 
is shown as a purple triangle. Note that no H2 gas velocity is 
shown in the lower panel due to its sparse detection along this 
slit. 



As for the other proposed hidden mass concentration, 
we do not see any prominent change in the stellar velocity 
dispersion or the molecular gas velocity along slit A or B 
(both roughly perpendicular to the original data in slit E) 
at the position proposed by TOO: the dispersion peak in the 
data of TOO, shown as slit E here, is therefore most likely a 
result of combined patchy extinction and young stars lead- 
ing to a variable stellar velocity dispersion depending on 
whether the bulge stars or young cluster stars dominate the 
light. 



6.2 A Nuclear Ring in M83 

Because there is no evidence to support a hidden mass con- 
centration at the centre of M83, we attribute the velocity 
gradients in the gas (ionised and molecular) to the presence 
of a shock; in fact, a shock and jump in the gas velocity on 
the bar's minor axis near the location of the north-east dust 
la ne is predicted by the n ucle ar ring and spiral si mulations 
of lRegan fc TeubenI (|2003 h and'M acieiewskif(|2004l ). Further- 
more, the simulations of iMacieiewskil predict shocks on the 
minor axis that are offset radially inwards to the maximum 
gas density, which we see is the case in Fig. |8] if extinction 
represents gas density. The second jump in the molecular 
gas velocity, 340pc along slit A, is also likely to be a shock: 
we can see in Fig. [7] that there is a significant quantity of 
dust 50pc to the south and south-west of the the nucleus, 
which may be spiralling inwards as in the nuclear spirals 
of Macieicwski (2004) • However, we estimate the posit ion 
of the kinematic centre defined in lSakamoto et al] (|2004l ) to 
be roughly 240pc along slit A which is almost equidistant 
from the two gradients in Vgas; therefore , given the scenario 
proposed bv lSofue fc Wakamatsul flOOi") a nd the s ymme try 
of the CO dust lanes and ring infSakamot o et al.l l|2004) . it 
is also possible for this gradient in Vgas to emanate from a 
shock on the inner edge of the south-western dust lane be- 
hind the bulge, rather than the aforementioned dust feature. 
We note that the H2 2-1 S(l) emission line can be excited 
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both by radiation (UV photons) and collisions (shocks) but 
it is most intense when coUisionally excited (compared to 
other lines). 

The CO bandhead, apart from being able to date the 
star clusters, also appears to be a useful diagnostic for esti- 
mating the timescale of the star formation for young clus- 
ters (< lOMyrs). This is due to the sharp onset of the fea- 
ture once the first red supergiants are formed: W[Ha] and 
W[PaQ] evolve relatively slowly in comparison and are not 
suitable for this purpose (Fig. [S] shows that the SSP models 
and the mixed model predictions for W[HQf] and W[Paa] are 
very similar). As discussed in i )3.3.2l we measure W[CO] to 
be consistently larger then the SSP model predictions based 
on the W[Paa] age. We found that the best way to rec- 
tify the conflict between W[CO] and W[Paa] was to replace 
the instantaneous burst with a finite episode of star forma- 
tion lasting 6 Myrs. This timescale is consistent with the 
collapse (and star formation) timescale for a giant molec- 
ular cloud with fiducial parameters in the simulations of 
iKrumholz et al.l (|2006l ). but it is also comparable to the dy- 
namical timescale for a half revolution around the galaxy 
centre at this radius. We also tested other finite episodes 
of star formation with varying star formation rates (expo- 
nential growth and exponential decay) and found our data 
couldn't reliably distinguish between them and the finite 
episode of constant formation. The real test of these differ- 
ent star formation scenarios is with clusters < 3Myrs old, 
where the model predictions difi'er most. We also tested a 
mixed population model of two SSP populations and found 
that two bursts separated by 5 Myrs with a mass ratio of 1:5 
almost reproduced the observed data trends in W[Paa] and 
W[CO], but slightly over predicted W[CO] for the youngest 
clusters. The observation of many more clusters in the north- 
east dust lane as an extension of the visible arc suggests that 
the stars are forming inside this dust lane and the mixing 
timescale of 6 Myrs reflects the formation timescale. How- 
ever, it is also feasible for the mixing to be a result of mul- 
tiple age populations. Although we were unable to create 
multiple burst models that fitted the data as well as the sin- 
gle burst lasting 6 Myrs, the scenario with star formation 
being triggered in the ring where the dust lanes join (on 
the minor axis of the bar) is appealing g iven that it has re- 
cently been fouiid in other nuclear ring s llAllard et al.ll2006l : 



.5 ^ 



ISarzi et al.ll2007l : iFalcon -Barroso ct al. '20071'). 

The CO maps of ISakam oto ct al. (20041) clearly show 
the molecular hydrogen in both dust lanes joining in a ring 
like structure. The material from the southern dust lane ap- 
pears to join the ring and move on the opposite side to the 
star burst arc to merge with the northern dustlane at the 
point where we see the youngest stars. Star formation at this 
point in the dust lane would usually be opposed because 
of the shear forces in action along the inevitable shocks; 
however, as i n the simulations of Athanassoula ( 1993) and 
iMacieiewskil l|2004 ). we observe the shock front (located by 
the gradient in Vgas) on the inner (southern) edge of the dust 
lane. On the outer edge, away from the shock front, where 
there is likely to be spurs and feathering (e.g. as found in 
iBonnell fc Dobbd l2006i l and less shear forces, star forma- 
tion may be less opposed. For stars formed here, on the 
western edge of the dust lane, orbital motion dictates their 
trajectory and they will propagate along X2 like orbits (il- 
lustrated in Fig. [TT] and also Fig. 15 of iRegan fc Teubeiil 
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Figure 11. A fiducial illustration of orbits (solid lines) and X2 
orbits (broken lines) in a ty pical bar pot ential (actually for the 
Milky Way), adapted from Binnev et al.l (:1991), with the false- 
colour image of Fig. [7] rotated and superimposed; the x-axis is 
aligned with the bar major axis while the y-axis is aligned with 
the bar minor axis. The bold line shows a 'cusped' xi orbit which 
separates the outer xi orbits from their smaller self-intersecting 
counterparts. Gas shocks along the leading edges of the bar and 
flows inwards to make the dust lane (at the point marked 1). This 
flow gradually funnels t he gas onto suitable X2 o rbits at the point 
marked 2, although see lReean fc TeubenI ll2003h and IMacieiewskil 
( 2004) for complete hydrodynamical simulations. If there was a 
source of star formation at this point (as proposed in the models 
of Regan & Teuben), such as the partner dust lane merging at 
the same location, then the age gradient along the circumnuclear 
arc could just be a result of the young stars moving on X2 orbits 
away from this location. 



l2003h to produce the 'star forming' arc. The gradient in age 
therefore refiects the older clusters having propagated fur- 
ther along their mutual X2 orbits. However, it is not obvious 
why the youngest star clusters we observe are 4.5 Myrs old; 
one would expec t them to be young er, unless the spraying 
of material as in iRegan et al.l (|l999l l led to star formation 
further up the dust lane. Equally, if star formation occurred 
on the opposite side to where we see the youngest clusters, 
triggered by the gas in the south-west of the ring merging 
with the south-west lane (the opposite of the previous sce- 
nario), then we would still expect an arc or ring of young 
stars to form but the timescales (dynamical and population 
synthesis) would agree better. If this were the case, we pre- 
sumably do not see the young stars in the north of the ring 
due to the low luminosity of very young clusters (< 5 Myrs) 
and the higher extinction there (seen in Fig. llOf) . Finally, 
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there may be star formation at both sites: the newly formed 
stars would still propagate on the same X2 orbits but the 
star formation from each point of contact does not need 
to be equal: the north-east d ustlane appears denser in the 
CO data of ISakamoto et al.l and could dominate the star 
formation. As previously discussed, we know that dual star 
formation events, separated by the dynamical timescale, can 
almost mimic a finite star formation episode of similar dura- 
tion, but our attempts to fit such a model to the data were 
not as successfuU as a single finite episode. 

A nuclear spiral or ring would cause the young clus- 
ters in the arc to be rotating around the galaxy centre in 
the same manner: clockwise as we see them in Fig. [2] (as- 
suming that the ring is co-planar with the galaxy). M83 
is almost face on, so the rotation is difficult to detect, but 
the slight inclination (24°. lTalbot et a"l]|l979l ) to our line-of- 
sight suggests that in the velocity data of TOO, we can see 
stars on the south-west side of the nucleus (coincident with 
the young clusters) streaming away from us and the stars 
on the north-east side stream towards us, in agreement with 
an inclined clockwise rotation. However, this slit (labelled 
slit E in Fig. [3)l mainly covers the background population 
rather than bright young clusters, so the velocity gradient 
could easily represent the overall rotation of the bulge. For 
consistency, let us estimate the dynamical time of the stars. 
From Fig. 1101 we see that the average velocity of the stars 
either side of the nucleus (and at the position of the arc) 
is roughly constant at ±25km s~^ in the rest frame of the 
galaxy nucleus along our line of sight. Assuming an inclina- 
tion of 24° then gives us an estimate for the circular velocity 
of 60km s"'^ ~ 60 pc Myr~^ (also comparable to the aver- 
age velocity dispersion). If the circumnuclear arc is around 
15" ~ 260pc in circumference, we expect the stars to take 
around 4 Myrs to travel from one end to the other. This is in 
good agreement with the age decrement found by compar- 
ing the cluster ages on either end of the luminous arc (~3.5 
Myrs). This point was also made in D06b, but the interpre- 
tation was significantly different. The dynamical timescale 
also changes with radius, which may help explain radial age 
gradients. 

We note that the northen dust lane does not appear 
to follow the same route as the stars along the arc in the 
HST images but instead continues to propagate south, past 
the nucleus. This might appear be a concern for the above 
scenarios. However, the CO maps (Sakamoto et al. 2004) 
clearly show the majority of molecular hydrogen gas follows 
the outer edge of th e circumnuc l ear ar c . We also no t e that 
in t he simulations of iPiner et al.l (119951 ) , iRegan et al.l (|l999l ') 
and iRegan fc TeubenI ('2OO3T similar features, where some 
gas appears to continue past the ring, can somet imes be seen, 
although they are very faint (not massive): in lRegan et al.l 
lll999h they are ref erred to as spray regions ( also discussed 
bv lSakamoto et al] ) and in lMacieiewskil (|2007l ) they are seen 
to be associated with regions of large shear. 



6.3 The displaced nucleus 

If the arc of young stars is a consequence of a nuclear spi- 
ral in the gas, we are left pondering why the only visible 
and kinematically detected nucleus at the centre of M83 is 
displaced from the centre of the galaxy and why we only 
see an arc of young stars rather than a ring. Apart from 



the arguments of ISofue~fc" Wakamatsul l| 19941 ) who suggest 
that a complete ring is inclined to the galaxy disk and ob- 
scured on the northern side, one should also consider a nat- 
ural m=l mode. Although rare, M83 is not the only spiral 
galaxy to show central asymmetries: NGC 1672 is barred spi- 
ral with a nucleus displaced from the kinematic and bulge 
centre, as D06b highlight; NGC 3504 appears to show a sin- 
gle dominant s piral arm in Ha and [Nil] encompassing an 
offset nucleus (|Emselleml l 200ll): NGC 3227 and NGC 7 130 
are both spirals highlighted in lMufioz Marfn et al] (|2007l ) as 
having off-centred nuclei; NGC 1808 is another barred spiral 
with a circumnuc lear starburst and nucleus offset fro m the 
kinematic centre l|Tacconi-Garman et al.l [19961 . |2005| ) while 
NGC 4303 shows asymmetric star formati on and dust lanes 
although the nucleus appears well centred l|Schinnerer et all 
l2002h . Although some of these galaxies are clearly interact- 
ing, none have clearly merged to an extent where the nuclei 
are likely to have been accre t ed ext ernally from the present 
interaction. iBournaud et al.l (|2005l ) have shown that m=l 
modes in the outer regions of disk galaxies can be repro- 
duced with a lopsided external accretion rather than a direct 
merger: if the outer regions of the galaxy are asymmetric, 
it may not be easy to define where the true morphological 
centre lies. That said, M83 does not look particularly dis- 
turbed or lopsided in the outer regions. If m=l modes in the 
outer disk can be induced by an external accretion event, 
it may also be possible for m=l modes to be induced by 
as ymmetric accret i on in the inner regions: in the CO maps 
of lSakamoto et al] (|2004l ). it appears that the gas inflow in 
M83 is slightly asymmetric, with the north-eastern dust lane 
denser than the south-west lane at the point of contact to 
the circumnuclear CO ring; this might cause unequal forc- 
ing on any nuclear spiral, although hydrodynamical models 
are need ed to confirm su ch a possibility. Finally, as noted 
by Sakamoto et alj l|2004h the nucleus could be suffering a 
natural, unprovok ed, instability as s een in the models of 
iTaga fc Ivd | 19981 ) or iLovelace et al.l (|l999l ) where the cen- 
tral nuclei are able to wander in the disk; such an instability 
might also be provoked by asymmetric inflow of gas. We 
noted earlier that the systematic velocity of the gas in the 
galaxy was taken to be 512 km s~^ but in Fig.[TD]the nucleus 
is moving towards us relative to this frame of reference at 
around 25km s~^, so not only is the nucleus displaced from 
the centre of the bulge, it is also decoupled kinematically. 

There is no evidence in the stellar kinematics of hidden 
mass concentrations at the locations previously proposed by 
TOO or D06(a,b). Without an integral field stellar velocity 
field of the central region of M83, we cannot completely rule 
out the presence of another hidden nucleus elsewhere, e.g. 
near the kinematical and bulge centre. However, given our 
extinction map (Fig. [J] from which we estimate Ay < 3 at 
the kinematic and bulge centre, corresponding to Ak< 0.3), 
the NICMOS F222M photometry (Fig. [2]), shows no massive 
luminous stellar feature other than the well known displaced 
nucleus: it is difficult to hide IO^Mq of luminous stars with 
only 0.3 magnitudes of extinction. 



7 CONCLUSION 

We have shown that the NIR indices W[PaQ] and W[CO] 
can be used to age-date young star clusters. Furthermore, 
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W[Paa] is much less affected by extinction than W[Ha], 
which requires significant correction. The young clusters in 
the circumnuclear arc are inconsistent with simple stellar 
populations produced by an instantaneous burst; the ob- 
servations require a finite episode of stax formation lasting 
around 6 Myrs. Such models together with our observations 
provide the clearest evidence yet of an age gradient along 
the star forming arc, as first found in HOI. However, our 
NIR observations further show that the circumnuclear arc 
does not stop at the edge of the dust lane, but continues 
through and along its north-west edge. 

Our stellar kinematics show no evidence for a second 
obscured nucleus at either the location of D06b or TOO and 
we conclude that the velocity gradients in the gas kinemat- 
ics (ionised and molecular) are a result of a shock front on 
the south-eastern edge of the north-eastern dust lane. We 
therefore explore alternative explanations for the formation 
of the circumnuclear arc, the most plausible being the X2 
orbital motion of the clusters from their place of birth in 
the dust lane. 

Future higher-resolution integral field spectroscopy in 
the NIR has the potential to build on this study and find 
the youngest clusters in the arc (suspected to be on the 
north-west edge of the northern dust lane) and also address 
previous claims of an age gradient perpendicular to the arc 
(HOI). 
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Table 4: The data for individual star forming knots. [1] The ID of the cluster in this paper. [2] Right ascention of the cluster (J2000). [3] Declination of 

the cluster (J2000) [4] The cluster ID in HOI. [5] The age of the cluster determined from the W[PaQ:] and W[CO]. [6] The age of the cluster from HOI. [7] 
The extinction Ay calculated from the HQ;:PaQ; decrement, given with 2cr errors (see §4.4). [8] Spatial contaminants: if the cluster is clearly resolved into 
multiple or blended components in HOI, we give the HOI ID and the age from that study in parentheses. [9] Spectral contamination: if the NIR velocity 
dispersion of the cluster is > 30km s~^, wc mark this column with a ^ . [10] Giant fraction: in the extraction of the kinematics of individual clusters, the 
best-fit kinematic template can be decomposed into a fraction of giant (III) stars and supergiant (I) stars; this is the fraction of giant stars in the best-fit 
template. [11] Supergiant fraction. [12] The Ha equivalent width of the cluster, uncorrected for extinction and with \a errors (but see §4.3). [13] The 
Pad equivalent width of the cluster, uncorrected for extinction and with \a errors. [14] The equivalent width of the CO bandhead, corrected for cluster 
dispersion and with la errors (but see §4.2). All: HOI cluster IDs are from their Table 2 unless accompanied with a superscript 3, indicating they are from 
their Table 3. 



